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ABSTRACT 

With thermal infrared observations detected by the NEOWISE project, we 
have measured diameters for 1742 minor planets that were also observed by the 
Infrared Astronomical Satellite (IRAS). We have compared the diameters and 
albedo derived by applying a spherical thermal model to the objects detected by 
NEOWISE and find that they are in generally good agreement with the IRAS 
values. We have shown that diameters computed from NEOWISE data are often 
less systematically biased than those found with IRAS. This demonstrates that 
the NEOWISE dataset can provide accurate physical parameters for the > 157,000 
minor planets that were detected by NEOWISE. 

1. Introduction 

Understanding the formation and evolution of the minor planets in our Solar System 
requires good knowledge of the physical parameters that describe each population of aster- 
oids and comets. Models of asteroid migration, for example, are critically dependent on 
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measurements of number counts, sizes and albedos. Visible light surveys have supplied most 
of the discoveries of asteroids and comets to date; however, these surveys tend to be biased 
against low albedo objects, and so they are underrepresented in number counts. Further- 
more, diameters derived from visible light observations are highly uncertain due to their 
linear dependence on albedo, which can vary by an order of magnitude for many types of 
asteroids. Radar observations, stellar occultations and in situ spacecraft imaging provide 
powerful means of obtaining precise diameter measurements, but these observations are lim- 
ited to a small subset of the ~500,000 minor planets known to exist today. While diameter 
and albedo can be measured more precisely than with visible light by applying radiometric 
models to objects for which thermal infrared data have been obtained, the size and albedo 
distributions obtained in this way will still be biased if the underlying source population was 
discovered by visible light surveys. In order to mitigate against these biases, it is necessary 
to undertake a survey that is capable of independent detection and discovery of asteroids at 
thermal infrared wavelengths. Yet the physical parameters such as size and albedo that are 
derived from such a new survey must be checked against other well-characterized datasets 
to ensure their reliability. 

The Wide-field Infrared Survey Explorer (WISE) is a NASA Medium-class Explorer mis- 
sion designed to survey the entire sky in fou r infrared wavelengths, 3.4, 4.6, 12 and 22 nm (de- 



noted Wl, W2, W3, and WA respectively) (IWright et al.ll2010l ; lLiu et al.ll2008l : Mainzer et al. 



20051 ). The final mission data products are a multi-epoch image atlas and source catalogs 



that will serve as an important legacy for future research. The portion of the pipeline dedi- 
cated to finding minor planets (called NEOWISE) has yielded observations of over 157,000 
minor planets, including Near-Earth Objects (NEOs), Main Belt Asteroids (MBAs), comets, 
Hildas, Trojans, Centaurs, and scattered disk objects ( jMainzer et al.ll2011al ). This represents 
an improvement of nearly two orders of magnitude more obje cts observed than W I SE's pre 



decessor mission, the I nfrared Astronomical Satellite (IRAS; iTedesco et al.lll988l ; iTedesco 



19921 ; iMatson. D.II1986I ). The WISE survey began on 14 January, 2010, and the mission 



exhausted its primary tank cryogen on 5 August, 2010. Exhaustion of the secondary tank 
and the start of the NEOWISE Post- Cryogenic Mission occurred on 1 October, 2010, and 
the survey ended on 31 January, 2011. 

In Mainzer et al. (hoilbh . we demonstrated that thermal models created for minor plan- 
ets using the WISE/NEOWISE dataset produce results that are in good agreement with di- 
ameters independently measured by radar, stellar occultations, and in situ spacecraft imag- 
ing. By taking the previously measured diameters from methods derived independently of 
infrared ther mal model assumptio ns, we were able to verify the accuracy of the color correc- 
tions given in IWright et al.l (120 10f ) for low temperatur e objects such as as teroids. We used 
the previously measured diameters combined with the IWright et al.l (120101 ) color corrections 
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in conjunction with a facete d spherical thermal model based on the Near-Earth Asteroid 
Thermal Model (jHarrislll998[ ) to predict the WISE magnitudes for ~50 asteroids; these were 
compared with the observed WISE magnitudes. We found that there were no systematic 
offsets between predicte d and observed mag nitudes for these objects, indicating that the 
color corrections given in IWright et al.l (120101 ) adequately describe the system response. For 
objects with WISE measurements in two or more bands with good signal-to-noise (for which 
the beaming parameter t] can be fit), we found that diameters can be determined to within 
±10%, and visible albedo pyto within ±20%. 

In this paper, we compare diameters and pv obtained with WISE/NEOWISE observa- 
tions for 1742 asteroids to those determined by IRAS and find that they are in generally good 



agreem ent with some slight systematic differences. This work, combined with iMainzer et al. 



( 1201 lbl ). demonstrates that the NEOWISE observations of minor planets will produce good 
physical parameters that will enable a wide range of scientific investigations. 



Observations 



We have assembled a list of objects that were observed by IRAS that NEOWISE de- 
tected dur ing the fully cryogenic portion of i t s mission. Of the ~2200 asteroids observed 



by IRAS (iTedesco et al.l Il988l : iTedesco I Il992l : Matson. D.lll986f ). we identified NEOWISE 



detections for 1742 objects. The observations of these objects were retrieved by querying 
the Minor Planet Center's (MPC) observation files to look for all instances of individual 
WISE detections of the desir ed objects that were r eported using the WISE Moving Object 
Processing System (WMOPS; IMainzer et al.ll2011al ). The resulting set of position/time pairs 
were used as the basis of a query of WISE source detections in individual exposures (also 
known as "Level lb" images) using the Infrared Science Archive (IRSA). In order to ensure 
that only observations of the desired moving object were returned from the query, the search 
radius was restricted to 0.3 arcsec from the position listed in the MPC observation file. Ad- 
ditionally, since WISE collected a single exposure every 11 seconds and observes each part of 
the sky an average of 10 times, the modified Julian date was required to be within 2 seconds 
of the time specified by the MPC. The following flag values were allowed: cc_flags = 0, P 
or p and ph_qual = A, B, or C. Objects brighter than W3 = 4 and W4 = magnitudes 
were assumed to have flux errors equivalent to 0.2 magnitudes due to changes to the shape 
of the point spread function as the objects became saturated, and a linear correction was 
applied to the W3 magnitudes in this brightness regime (the WISE Explanatory Supplement 
contains a more detailed explanation). Per the Explanatory Supplement, objects brighter 
than W3 = —2 and W4 = —6 were not used. Each object had to be observed a minimum 
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of three times in at least one WISE band, and it had to be detected at least 40% of the 
time when compared to the band with the maximum number of detections (usually, though 
not always, W3). The WMOPS system is designed to reject inertially fixed objects such as 
stars and galaxies in bands W3 and W4. Nonetheless, the individual images at all wave- 
lengths were compared with WISE atlas coadd and daily coadd source lists to ensure that 
inertially-fixed sources such as stars and galaxies were not coincident with the moving object 
detections. This check is particularly important in bands Wl and W2 where the density of 
background objects (and hence the probability of a blended source) is higher than at longer 
wavelengths. Any remaining blended sources in bands Wl and W2 were removed. Some 
objects were observed at multiple epochs, and observations separated by more than three 
days were modeled separately. 



3. Thermal Model and Reflected Sunlight Fits 



In the Standard Thermal Model (STM) of Lebofsky fc Spencer (jl989l ). the temperature 
of an asteroid is assumed to be maximum at the subsolar point and zero on the point 
opposite of this; this is the ca se of an object with zero thermal inertia. In contrast, in the 
Fast Rotating Model (FRM) (ILebofskv et al.lll978[ IVeeder et al1ll989l : iLebofskv fc Spencer 
19891 ). the asteroid is assumed to be rotating much faster than its cooling time, resulting in a 
constant surface te mperatur e acros s all longitudes. The so-called "beaming parameter" (rj) 
was introduced by iLebofskyl (119861 ) in the STM to account for the enhancement of thermal 
ra diation obser ved at small phase angles. The near-Earth asteroid thermal model (NEATM) 
of lHarrisI (119981 ) also uses the beaming parameter r] to account for cases intermediate between 
the STM and FRM models. In the STM, 7] is set to 0.756 to match the occultation diameters 
of (1) Ceres and (2) Pallas, while in the FRM, t] is equal to tt. With NEATM, 77 is a free 
parameter that can be fit when two or more infrared bands are available. 

We modeled each object as a set of triangular fac ets covering a spherical surface with 
diameter equal to the ground-truth measureme nt fc.f. iKaasalainen et al.l 120041 ) . The tem- 
perature for each facet was computed, and the Wright et al. ( 201ol ) color corrections were 
applied to each facet. The emitted thermal flux for ea ch facet was c a lculat ed using NEATM 
along with the bandcenters and zero points given in IWright et al.l (120101 ); nightside facets 
were assumed to contribute no flu x. The emissivity, e, was assumed to be 0.9 for all 
wave lengths (c.f. 



Harris et al. 



20091) . The objects' absolute magnitudes (H) were taken 



from lWarner. Harris &: Praved ( 120091 ) when available; otherwise, the values were taken from 
the MPC's orbital element files. Unless a direct measurement of G was available from 



Warner. Harris &: Praved (120091 ) . we assumed a G value of 0.15. In general, minor planets 
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detected by NEOWISE in bands Wl and W2 contain a mix of reflected sunlight and ther- 
mal emission. Thus, it was necessary to incorporate an estimate of reflected sunlight into 
the thermal model in order to use data from bands Wl and W2. In order to compute the 
fraction of reflected sunlight in bands Wl and W2, it was also necessary to compute the 
ratio of the albedo at these wavelengths compared to the visible abl edo (pm/py); t he sola r 
spectrum was approximated as a 5778 K blackbody. As described in lMainzer et al.l ( 1201 lbl ). 
we assumed that pm = pzAy,m — Vi.&nm- The flux from ref lected sunlight was computed 
for each WISE band using the IAU phase curve correction of iBowell et al.l (119891 ). Thermal 
models were computed for each object by grouping together observations no more than a 
three day gap between them. 

Error bars on the model magnitudes and subsolar temperatures were determined for 
each object by running 25 Monte Carlo (MC) trials that varied the objects' H values by 
0.3 magnitudes and the WISE magnitudes by their error bars using Gaussian probability 
distributions. The minimum magnitude error for all WISE measurements fainter than W3 = 
4 and W4 = 3 mag nitudes was 0.03 magnitudes, per the in-band repeatability measured in 
Wright et al.l ( 120101 ). For objects brighter than W3 = 4 and W4 = 3, the error bars were 



increased to 0.2 magnitudes, as these magnitudes represent the onset of saturation (see the 
WISE Explanatory Supplement). For those objects for which r\ and pm/pv could not be 
fitted, rj was set to 1.0 and allowed to vary in the MC trials by 0.25, and pm/pv was set to 
1.4 and allowed to vary by 0.5. These default values were used because they are the mean 
and standard deviation of rj and Pm/pv for the objects for which these parameters could be 
fit. The error bar for each object's model magnitude was equal to the standard deviation of 
all the MC trial values. 

Although many asteroids are known to be non-spherical, the WISE observations gen- 
erally consisted of ~10-12 observations per object uniformly distributed over ~36 hours 
( jWright et eil.l 12010c iMainzer et al.l 1201 lal ). so on average, a wide range of rotational phases 
was sampled. Although the variation in effective spherical diameter resulting from rotational 
effects tends to be averaged out, caution must be exercised when interpreting effective di- 
ameter results using spherical models for objects that are known to have large-amplitude 
lightcurve variations. In our sample, 179 objects have NEOWISE observations at multiple 
epochs (meaning that the groups of observations were separated by more than 10 days). 
Of these, all but 24 had diameters that agreed to within 10%, and most of the remaining 
objects had W3 peak-to-peak amplitudes >0.3 mag, indicating that they are likely to be 
non-spherical. Non-spherical objects observed at different viewing geometries can lead to 
different diameters when using a spherical model. However, since 87% of the objects with 
multi-epoch observations have diameters that agree to within 10%, we conclude that multi- 
ple observational epochs do not contribute significantly to the differences observed between 
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IRAS and NEOWISE diameters. 



Results 



Tedesco et al.l (120021 ) reported observations of ~2200 asteroid within the Supplemental 



IRAS Minor Planet Survey (SIMPS). The diameters and albe dos in that wor k were computed 
using the Standard Thermal Model. Similar to the work of IWalkerl (120031). who computed 
NEAT M diameters and albedos for 654 asteroids using IRAS fluxes, iRyan &: Woodward 
( 120101 ) applied a spherical NEATM model to ^1500 aster oids found in the SIMPS and 
Mid-Course Space Experiment (MSX; iTedesco et al.l l2004bl ) data. Figures [T] and H] show 



the compariso n betwe en di ameters found using NEOW ISE observations and those given by 
Tedesco et all J2004ah and iRvan & Woodward! (boioh . respectively. Figures [3] and @] show 
the comparison betwe en NEOWISE and Tedesco/Ryan albedos. The diameters given by 
Tedesco et al.l (j2004al ) are systematically lower t han the NEOWISE diameters, while the 
dia meters computed by IRvan fc Woodward! (120 !0l ) are systematically higher. Also, py given 
by Tedesco et al.l (2004a ) is systematically higher than the NEOWISE py values, but the 
Ryan fc Woodward! (120101) value s te nd to be slightly lower. In Figure we show the di- 
ameters of ITedesco et al.l ( 12004a! ) and lRyan fc Woodward! (120101 ) compared to the diameters 
of 78 objects with di ameters measured by radar, occultation, or spacecraft imaging (see 
Mainzer et al.l (12011b! ) for a complete list of references). These objects have been selected 
to have W3 peak-to-peak amplitudes <0.5 magnitudes to avoid the worst complications 
caused by applying spher i cal mo dels to highly elongated objects. As can be seen in the 
figure, the ITedesco et al.l (|2004al ) diameter s are skewed toward sligh tly smaller sizes than 
the radar, etc. measurements, whereas the IRyan fc Woodward! (120101 ) diameters are biased 
toward larger sizes. We propose that these biases are caus e d by the band-to-band color cor- 
rections derived for IRAS. As described in ITedesco et al.l (120021 ). band-to-band corrections 
were derived by requiring that the 10 and 20 /zm IRAS observations simultaneously matched 
a diameter derived from a stellar occultation of (1) Ceres. Because the STM parameters 



were set to match Ceres, ITedesco et al.l ( 120021 ) found a 7% difference between diameters de- 
rived for 1 3 objects with diame t ers m easured from stellar occultations. The NEATM model 
applied by IRvan fc Woodward! ( 120101 ) does not appear to address this issue. Furthermore, 
Rvan fc Woodward! J201(F applied cutoffs to 77 of 0.75 and 2.75. As described above, in our 
NEATM models applied to the NEOWISE data, we cut off r\ at the limit predicted by the 
Fast Rotating Model, it. In the STM, rj was set to 0.756 by matching the diameters of Ceres 
and Pallas; however, there is no reason in principle why it cannot range to lower values. In 
our application of NEATM, we allow r] to range arbitrarily low and find no fitted 7] values 
less than ~0.53. 
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As described above, in lMainzer et al.l (l2011bf ). we used asteroids with previously known 
diameters and H magnitudes to compare our model and observed WISE magnitudes, and 
we found no systematic offsets. We conclude that the diameters for minor planets derived 
from NEOWISE are in generally good agreement with those found by IRAS and are lik ely 
more free of systematic bi ases than the di ameters provided in eit her iTedesco et al.l (120021 ) or 



Ryan fc Woodward! ( 120101 ) . Together with lMainzer et al.l (j2011bl ). this demonstrates that the 



NEOWISE dataset will produce good quality physical parameters for the > 157,000 minor 
planets it contains. 
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Fig. 1. — : Comparison of the difference in diameters given by iTedesco et al.l (J2004aj) and 
the diameters found by applying a spherical NEAT M model to NEOWIS E measurements for 
1742 asteroids. It can be seen that the values from lTedesco et al.l (|2004al ) tend to be slightly 
smaller than NEOWISE-derived diameters. The NEOWISE diameters diverge widely from 
IRAS for the smallest objects; however, these objects are observed by NEOWISE with high 
signal-to-noise ratio. 
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Fig. 2. — : Comparison of the difference in diameters given by lRyan fc Woodward! (120101 ) and 
the diameters found by applying a spheri cal NEATM model to NEO WISE measurements for 
1155 asteroids. The diameters found by iRyan fc Woodward! (120 lOl ) tend to be larger than 
those found by NEOWISE. 
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Fig. 3. — : Comparison of the difference in albedos given by lTedesco et al.l (j2004al ) and the 
diameters found by applying a spherical NEAT M model to NEOWISE measurements for 
1742 asteroids. It can be seen that the values from lTedesco et al.l (]2004al ) tend to be slightly 
higher than NEOWISE-derived albedos. 
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Fig. 4. — : Comparison of the difference in albedos given by iRyan fc Woodward! ( 120101 ) and 
the albedos found by applying a sph erical NEATM model to N EOWISE measurements for 
1155 asteroids. The albedos found by lRyan fc Woodward! (120 lOl ) tend to be lower than those 
found by NEOWISE. 
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Fig. 5. — : Comparison of diameters found by radar, stellar o ccupation, or i n situ s pace- 
craft imaging for 78 ast eroids and the diameters found by iTedesco et all (j2004al ) and 
Ryan &: Woodward! (120101 ) . The objects were selected to have W3 peak-to-peak amplitudes 
<0.5. It can be seen that both methods of computing IRAS diam eters are somewhat biased 
when compared to the radar diameters. In lMainzer et all ( 1201 lbl ). we show that there is no 
systematic bias found in model and observed WISE magnitudes for a similar set of objects 
with independently measured diameters. 



